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Abstract 

Wigner molecules formed at high magnetic fields in circular and elliptic quan- 
tum dots are studied by exact diagonalization (ED) and unrestricted Hartree-Fock 
(UHF) methods with multicenter basis of displaced lowest Landau level wave func- 
tions. The broken symmetry states with semi-classical charge density constructed 
from superpositions of the ED solutions are compared to the UHF results. UHF 
overlooks the dependence of the few-electron wave function on the actual relative 
positions of electrons localized in different charge puddles and partially compen- 
sates for this neglect by an exaggerated separation of charge islands which are more 
strongly localized than in the exact broken-symmetry states. 
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At high magnetic fields the electron systems in circular quantum dots form 
Wigner molecules [1,2] in the internal structure of the system. Deformation 
[3,4] of the circular symmetry allows the molecules to appear in the labora- 
tory frame only [4] when the classical counterpart [5] of the quantum sys- 
tem possesses a single lowest-energy configuration. Otherwise the charge den- 
sity at high field is a superposition of equivalent semi-classical densities and 
the quantum system undergoes symmetry transformations when the magnetic 
field is increased [4]. These transformations are associated with level cross- 
ings at which the ground state is two-fold degenerate. Superposition of the 
states of the degenerate levels allows [4] to extract the semi-classical broken- 
symmetry charge density into the laboratory frame. On the other hand the 
unrestricted Hartree-Fock (UHF) produces [1] broken-symmetry states for 
Wigner molecules. In the infinite magnetic field limit UHF gives exact [6] 
results for the total energy. At finite magnetic fields for which exact broken- 
symmetry eigenstates exist the artifactal symmetry breaking cannot be blamed 
to the inaccuracy of the UHF [4]. In this paper we look for the effects neglected 
by UHF comparing the ED and UHF solutions for elliptical and circular dots. 

We assume a spin-polarization of electrons at high magnetic field (0, 0, B) 
oriented perpendicular to the quantum dot plane and use the Landau gauge. In 

the ED, described in detail in Ref. [4], the single electron wave functions used 
for construction of the Slater determinants are obtained via diagonalization of 
the single electron Hamiltonian in the multicenter basis [7,8,9] of M displaced 
lowest Landau level wave functions 



where a is treated as variational parameter. In the present UHF approach the 
one-electron orbitals (1) are optimized self-consistently. We study up to A^" = 4 
electrons, use the material data of GaAs and a basis of 12 centers {xk, Vk) put 
on an ellipse with a size determined variationally. Basis (1) of displaced lowest 
Landau level wave functions reproduces [4] also higher Fock-Darwin bands. 
Contrary to previous multicenter HF calculations [6,9] using a single center 
per electron, the present HF approach produces results which are exact in the 
UHF [10] sense. 

Classical system of three electrons in an elliptical confinement potential with 
hujx = 3 meV and huty = 4 meV possesses two equivalent lowest-energy config- 
urations [cf. inset of Fig. 1] and the quantum system undergoes parity trans- 
formations [4] with the magnetic field [cf . Fig. 1] . Superposition [4] of the two 
lowest-energy eigenstates 
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yields a brokcn-symmctry (BS) charge density with a distinct electron separa- 
tion. Fig. 1 shows that in contrast to the exact ground-state energy the UHF 
energy estimate is a smooth function of the magnetic field. 

The charge densities of considered states are shown in Figs. 2(a) and 2(b) 
for two magnetic field values corresponding to the even-odd energy crossing 
presented in Fig. 1. The phase in the BS state is chosen such that the 
electrons are localized at the classical Wigner molecule positions. Notice that 
in the UHF the separation of electrons is more pronounced than in the exact 
BS states. Fig. 2(c) shows the pair correlation function (PCF) [2] for the UHF 
and the exact BS state corresponding to the charge density of Fig. 2(a) with 
the position of one of the electrons fixed at two different locations: in the 
center and on the edge of the central charge puddle. In contrast to the exact 
BS state in the UHF wave function the two electrons are insensitive to the 
actual position of the third electron in its charge puddle. This is a consequence 
of the single-determinantal form of the UHF wave function, and can be easily 
explained for two electrons. In the spin-polarized two electron Wigner molecule 
the UHF spatial wave function is given by ^'a(ri)^'/3(r2) — ^'y3(ri)^'a(r2), where 
l^'ap and I^^P are the charge densities of separate charge puddles a and 
/3. Wave functions and are orthogonal due to the vanishing overlap 
between the puddles. The calculation of the PCF gives (up to a constant) 
PCF(r„,r5) = |^«(r„)|2|^'^(rfe)|2 + \^ fi{va)\^ . For rj, inside puddle /3 
the second term of the sum vanishes and the remaining one signifies that the 
probability of finding an electron in point inside puddle a is independent 
of the actual position of the second electron in puddle /3. 

The UHF self-consistency is reached only in one of the two classical orienta- 
tions [cf. inset of Fig. 1] in which the UHF energy is minimal. On the other 
hand, the exact BS states can be oriented under an arbitrary angle [cf. Fig. 
2(d)] by modifying the phase in Eq. (2). Moreover, since the BS state is con- 
structed with states of opposite parities all the plots in Fig. 2(d) correspond 
to the same value of the kinetic, potential and electron-electron interaction 
energies equal to the arithmetic average of the expectation values for '^odd 
and '^even states. 

Fig. 3 shows the two lowest energy levels and the UHF total energy calculated 
with respect to the lowest Landau level for the elliptical dot with h,LOx = 3 
meV and hujy = 4 meV. For these values the classical counterpart of the four- 
electron system is unique and conform with the symmetry of the confinement 
potential, so that the Wigner crystallization is visible in the exact quantum 
ground-state for an arbitrary magnetic field after the MDD decay. In this case 
the MDD decay is a continuous process and appears at the anticrossing around 
6 T. 

The inset to Fig. 3 shows the charge densities calculated with the UHF and 
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ED methods. In ED the charge density in between the charge maxima takes 
on larger values than in UHF in which the separation of electrons is more 
distinct like ior N = 3 [cf. Fig. 2(a) and (b)]. 

The average radius of the charge puddle as obtained for the two-electron 
ground state in the circular quantum dot {huJx = = 3 mcV) is displayed 
in Fig. 4. The ED and UHF values are similar below the MDD breakdown 
{B < 5.6 T). The exact value present discontinuities at the angular momen- 
tum transitions. After the MDD decay the UHF value is close to the average 
around which ED and BS results oscillate, but at higher fields it becomes an 
upper bound for these oscillations. The inset of Fig. 4 shows the exact BS 
and UHF charge density for B = 11.8 and 28.9 T. BS charge densities for 
the two values of the magnetic field have been obtained from superpositions 
of the degenerate states with angular momenta -5,-7 and -9,-11 (in h units), 
respectively. The charge maxima in the exact BS solutions are less strongly 
separated. For B = 28.9 T the charge density islands of the exact BS have 
crescent shape while the UHF charge puddles are more oval and the distance 
between them is larger. 

The difference in shape of the separated charge islands in circular dots is 
largest ior N — 2. For larger the charge puddles in the exact solutions 
are less spectacularly spread. Fig. 5 shows the comparison of the exact BS 
and UHF charge densities for = 3. The BS plots for A^ = 3 correspond to 
degeneracy of states with angular momentum —9h, —12% for B = 7.5 T and 
-12h,-15h for 5 = 15 T. 

In summary, we have presented a comparison of the UHF and the ED results 
for the charge density of Wigncr molecules in circular and elliptical dots. For 
the comparison we have used the broken-symmetry states obtained from the 
superposition of the exact eigenstates. We have found that the UHF exagger- 
ates the separation of the electron charge densities in the laboratory frame. 
In this way the UHF method partially compensates for the overlooked cor- 
relations related to the reaction of electrons on their actual position in the 
separated charge density islands. This reaction is of smaller importance for 
larger magnetic fields, for which the charge density islands shrink to points, 
which explains the vanishing of the UHF energy overestimation in the infi- 
nite magnetic field [6]. Due to the exaggerated electron separation the charge 
islands forming Wigncr molecules calculated in the UHF shrink with the mag- 
netic field faster than in the exact broken-symmetry states. 
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Fig. 1. Energy of the lowest even (dashed curve) and odd parity (sohd curve) levels, 

the broken symmetry state (dotted curve) and UHF energy calculated for = 3, 
HlOx = 3 meV and fiLOy = 4 meV. Inset shows the two equivalent classical configu- 
rations. 



5 



B=8.75 T 
even ^ odd 5 




BS - = UHF 




B=17.25T 

7t/8 , 71/4 - 




Fig. 2. (a) The charge densities for the even, the odd, the broken-symmetry (BS) 
state and the UHF for N = 3, hco^c = ^ meV, tiLOy = 4 meV at S = 8.75 T. (b) Same 
as (a) at B = 17.25 T. (c) Pair correlation function plots for UHF and BS state 
for an electron fixed at (0.13.4 nm) [left panel] and (-3 nm,10.4 nm) [right panel] 
marked by crosses, (d) BS states for given shifts of phase (j) [Eq. (2)] with respect 
to the BS plot in (b). Plots (a), (b) and (d) have the same scale given by the length 
bar in (a). 
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Fig. 3. Lowest energy levels for A'' = 4, Tiuj^ = 3 meV and hujy = 4 meV (solid 
lines) and total UHF energy (dotted) line. Inset: Charge density calculated in the 
UHF (two upper plots) and the ED (two lower plots) for S = 8 T (left side) and 
5 = 20 T (right side). 
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Fig. 4. Average radius of the charge puddle (< >^/^) for two electrons in a circu- 
lar quantum dot {hiOx = htOy = 3 meV). Dotted, dashed and solid lines correspond 
to UHF, BS, and ED results. Inset: charge density obtained in the UHF (two upper 
plots) and in the exact broken-symmetry state (two lower plots) for B=11.8 T (left 
side) and 28.9 T (right side). 
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Fig. 5. Charge density for = 3 electrons in a circular quantum dot 
{hujx = fioJy = 3 meV) calculated with UHF (upper plots) and for exact BS states 
(lower plots) for 5 = 7.5 T (left plots) and B = 15 T (right plots). 
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